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Введение. Статья посвящена теоретическому описанию 
процесса переноса ионов в мембранных системах в галь- 
ваностатическом режиме. В качестве мембранной системы 
рассматривается канал обессоливания электродиализного 
аппарата. Цели работы: создание и верификация двумер- 
ной математической модели стационарного переноса 
ионов соли в канале обессоливания электродиализного 
аппарата для гальваностатического режима. 

Материалы и методы. Предложена новая модель перено- 
са ионов. Она основана на системе уравнений Нернста — 
Планка — Пуассона для электрического потенциала и на 
уравнении для функции электрического тока. Получено 
численное решение краевой задачи модели методом ко- 
нечных элементов с помощью программного пакета 
Сот50о1 МиШрйу51с$. 

Результаты исследования. Разработанная математическая 
модель позволяет описать стационарный перенос ионов 
бинарной соли в канале обессоливания электродиализного 
аппарата. При этом учитываются нарушение электроней- 
тральности раствора и формирование расширенной обла- 
сти пространственного заряда при сверхпредельных токах 
в гальваностатическом режиме. Об адекватности постро- 
енной модели свидетельствует хорошее совпадение физи- 
ко-химических характеристик переноса, рассчитанных по 
моделям для гальваностатического и потенциостатическо- 
го режимов. 

Обсуждение и заключения. Разработанная модель позво- 
ляет интерпретировать результаты экспериментальных 
исследований переноса ионов в мембранных системах, 
если данный процесс протекает в гальваностатическом 
режиме. Некоторые электрокинетические процессы связа- 
ны с появлением расширенной области пространственного 
заряда при сверхпредельных токах. Описывая формирова- 
ние указанной области, можно выяснить, каким образом 
зависящие от нее процессы влияют на перенос ионов в 
гальваностатическом режиме. 
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Введение. Мембранные системы составляют основу электродиализных аппаратов, нано- и микрофлю- 
идных устройств, которые применяются при очистке воды, переработке сельскохозяйственной продукции (мо- 
лока, вина и т. д.), выполнении химических анализов и в других сферах деятельности [1-4]. В многочисленных 
математических моделях процессов массопереноса в мембранных системах для потенциостатического или по- 
тенциодинамического режимов электрический режим задается как скачок потенциала между двумя эквипотен- 
циальными плоскостями, параллельными мембранам. Подробный обзор работ, посвященных моделированию 
для потенциостатического режима, представлен в [5—7]. 

В то же время в практике электродиализа, электрохимической характеризации мембран (хронопотен- 
циометрия, импедансометрия и др.) часто используется гальваностатический режим, при котором на межфаз- 
ной границе поддерживается постоянная средняя плотность тока. Об этом режиме собран огромный объем экс- 
периментальных данных, которые необходимо интерпретировать [8—10]Ошибка! Закладка не определена.. 
Исследования в области математического моделирования гальваностатического режима ведутся по нескольким 
направлениям. 

Первое направление — метод обратной задачи. Как ясно из названия, речь идет о решении обратной 
задачи: для заданной плотности тока на межфазной границе «раствор — мембрана» находится соответствую- 
щий скачок потенциала, а далее рассматривается задача для потенциостатического режима [11]. Низкая эффек- 
тивность данного метода обусловлена тем, что его реализация требует многократного решения задачи в потен- 
циостатическом режиме для одного заданного значения плотности тока. 

Второе направление — метод декомпозиции. При этом система уравнений Нернста — Планка — Пуас- 
сона заменяется системой декомпозиционных уравнений [12—16]. Предположение о квазиравномерном распре- 
делении заряда позволяет получить модель для гальваностатического режима в приближении закона Ома [17— 
20]. 

Третий подход можно назвать прямым методом. В этом случае для плотности тока в канале обессоли- 
вания выводится уравнение, заменяющее уравнение Пуассона [21]. 

Гальваностатический режим можно описать иначе — с помощью численного решения уравнений 
Нернста — Планка — Пуассона для электрического потенциала со специальным граничным условием, которое 
позволяет установить плотность тока как параметр, задающий электрический режим в системе. В [22, 23] для 
одномерного случая производная по времени градиента электрического потенциала определялась как явная 
функция плотности тока. Это отличает подходы авторов указанных работ от потенциостатических моделей, в 
которых задается разность потенциалов. 

В данной статье представлена стационарная модель процесса переноса ионов в мембранных системах 
для гальваностатического режима. Она основана на системе уравнений Нернста — Планка — Пуассона с гра- 
ничным условием, которое позволяет установить плотность тока как параметр, задающий электрический режим 
в системе. В этом решение аналогично [22, 23]. Отличие же состоит в том, что предлагаемая модель двумерна и 
учитывает непостоянство плотности тока по длине канала. 

Материалы и методы. Под мембранной системой подразумевается канал обессоливания электродиа- 
лизного аппарата (ЭДА), образованный анионообменной (АОМ) и катионообменной (КОМ) мембранами. Через 


него со средней скоростью У прокачивается раствор бинарного электролита. На рис. | х — нормальная к по- 


верхности мембраны координата, изменяющаяся от 0 (граница с АОМ) до й (граница с КОМ); у — тангенци- 
альная к поверхности мембраны координата, изменяющаяся от 0 (вход в канал) до [ (выход из канала). 
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Рис. 1. Схема канала обессоливания ЭДА. Показаны концентрационные профили катионов С\ (сплошная линия) и 


анионов С› (пунктирная линия). Скорость вынужденного течения У 


Система уравнении. Рассмотрим двумерный стационарный случай системы уравнений, описывающей 
перенос бинарного электролита в отсутствие химических реакций [24]: 


Е. Е а 

Л 51. С; УФф-Р; УС; + СУ, У ВАА (1) 
КТ 

—@У 7 =0, 1=Ь2, (2) 
508,Аф =-Р(аС1 + 52С.), (3) 
ТЕР +22),). (4) 
Здесь |, О; зи С, — соответственно поток, коэффициент диффузии, зарядовое число и молярная концен- 
трация {1-го иона; ф — электрический потенциал; У — скорость потока раствора электролита; &0 — электриче- 


ская постоянная, &, —— относительная диэлектрическая проницаемость раствора электролита (предполагается 


постоянной); [ — плотность тока; ЕЁ — постоянная Фарадея; К — газовая постоянная; ТГ — абсолютная темпе- 


ратура; Л, /2, 1, Ф,С., С› — неизвестные функции координат х иу. 


Уравнения Нернста — Планка (1) описывают поток ионов, обусловленный миграцией в электрическом 
поле, диффузией и конвекцией; (2) — уравнение материального баланса в стационарном случае; (3) — уравне- 
ние Пуассона для потенциала электрического поля; (4) — плотность тока в растворе электролита. Будем счи- 
тать, что распределение скорости в канале соответствует течению Пуазейля [24]: 


у. =0, У, =6% > а | (5) 


Для определения перечисленных выше неизвестных функций необходимо задать краевые условия для 
системы (1)—3). Рассмотрим два электрических режима: потенциостатический, когда скачок потенциала на 
границах системы задается постоянным, и гальваностатический, когда постоянна плотность тока, протекающе- 
го через межфазную границу. 

Граничные условия для моделирования потенциостатического режима. Будем считать поверхности 
ионообменных мембран эквипотенциальными. Система (1)—(4) включает потенциал электрического поля толь- 
ко в форме производных по пространственной координате. В этом случае существенным является только ска- 


чок потенциала Аф = ф<(й, у) — Ф<(0, у) ‚ где Аф — известная функция, поэтому положим, например: 
<(0, у)=0. (6) 
Тогда 
До = соп5Е . (7) 
Условия (6) и (7) определяют потенциостатический режим. Остальные граничные условия приводятся 
ниже. 


На межфазной границе «АОМ — раствор» (х = 0) концентрация коионов (катионов) определяется из 
условия непрерывности их потока у границы «мембрана — раствор» с учетом селективных свойств АОМ [16]: 
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[7 Е =\ о (1-7>)1, (0,5) (8) 


х ВТ ах эВ 


где Т; (= 1, 2) — эффективные числа переноса противоионов в мембране (КОМ и АОМ соответственно); Т; — 


числа, близкие к 1, причем для идеально селективной мембраны /; = 1, и условие (8) превращается в условие 


непроницаемости мембраны для коионов. 
Концентрация противоионов (анионов) зависит от обменной емкости КОМ, что можно задать в виде: 
С. (0,у) =С.„ =М,С.. (9) 
Здесь постоянная №, показывает, во сколько раз эта концентрация отличается от концентрации в объеме раство- 
ра [25]. 
На межфазной границе «раствор — КОМ» (х = й) для концентраций ионов приняты условия, аналогич- 
ные условиям на границе «АОМ — раствор» (х = 0): 


С(, у) = Ст = № Со, (10) 
ОС Е д 1-7) Г. (®, 
+ м (у) = (27215 >) :) х( У (10) 
Ох ЮТ Ох <о ЕР 
На входе в канал (у = 0) принято равномерное распределение концентрации ионов: 
С: (х.0) = Су, #=12. (12) 


Условие для потенциала электрического поля получено из уравнений (1) и (4) с учетом отсутствия про- 
текания тока через вход /, (х, 0, В) =0: 


Оф(х,О КТ ОС1 (х,0 ОСо (х,О 
оо Е ЗЕ [а РН, а] (13) 
ду ЕГО + &22.5)Со ду ду 
На выходе из канала (у = Г) ионы свободно выносятся потоком раствора: 
. = 5 2 - | 
(п, ]; ) = = 50; С, Уф-Б; УС; +СМ) = СИ, 1=Ь2. (14) 


Система (14) означает также, что сумма диффузионной и миграционной тангенциальных составляю- 
щих потоков анионов и катионов равна 0: 


бе ЙО 
| _ итас, |) 0,5=1,2. (15) 


Для потенциала принято «мягкое» условие, означающее отсутствие резких изменений потенциала на 





выходе из канала: 


= 0. (16) 
ду 


Краевая задача, включающая уравнения (1)—(4) и краевые условия (6)—(16), моделирует потенциостати- 
ческий режим, причем его определяющее ключевое условие — (7). 

Граничные условия для моделирования гальваностатического режима. При моделировании гальвано- 
статического режима условие (6) остается, а (7) необходимо заменить условием, связанным с заданным значе- 


нием средней плотности тока 1, на межфазной границе «раствор — КОМ» (х = Й). Для вывода такого гра- 


ничного условия подставим соотношения (1) в (4) и выразим градиент потенциала электрического поля: 
АТ = - 
Уф = — р р. ([ + Е(аПУС я 520>У С) == Е(аС = 52С5)У). (17) 
Е (3 ВС -Ё 5210№С.) 
Полагая в (17) Хх = ‚, получаем соотношение, связывающее градиент электрического потенциала с за- 


данным значением плотности тока на границе, то есть граничное условие на межфазной границе «раствор — 
КОМ»: 


ас, вс, 
Г; + РО — + Е 5205 Е 
Ви = ВКТ дх ах 


(й, у) (18) 
Ох Е2 2 О:С! ЕН &>02С. 


При этом плотность тока /,. должна удовлетворять условию: 
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1! | 
- [1.4 = 4. (1) 
0 


Математическая модель гальваностатического режима состоит из системы уравнений (1)—(4). Гранич- 
ные условия (18) и (19) заменяют условие (7). Остальные граничные условия совпадают с условиями для по- 
тенциостатического режима. 

Скачок потенциала в гальваностатическом режиме является вычисляемой величиной. 

Преобразование граничных условий для моделирования гальваностатического режима. Условие (19) 
неудобно для численного решения, так как содержит интеграл. Ниже приводится один из вариантов преобразо- 


вания этого условия. 


В стационарном случае плотность тока 1 является соленоидальным вектором. Действительно, если 


умножить (2) на 5; и сложить, то ЧУ/ = 0. Следовательно существует такая функция 1, что 


Е, 


| 20 
9х ” ду ь — 
С использованием функции 1 условия (18) и (19) переписываются следующим образом: 
| ы = + Ес О, ет + Ре,Ь, — 
х х 
ААУ) ===> 2 2 (п, у), (20) 
Ох Е Са. С, 
1! 1-0 1 | 
= [1 =-- [уу = @)-п@0)) =. (22) 
в. Г. ду | 
Уравнение (22) можно переписать: 
(й,1)-п(ф,0) ==. (23) 


Чтобы замкнуть систему выражений, необходимо получить уравнение для функции 1. С этой целью, 


как и в работах [15, 16], введем в рассмотрение линейный дифференциальный оператор, который является 
функцией вихря (ротором) в двумерном случае, для произвольного двумерного вектора У: 
и _ СУ 


И) = —* | 
(и) х в (24) 





Несложно проверить, что: 
1) г(Уи) = 0 для любой гладкой функции И; 
2.) (и) = (Уи,\)1 + иг(’) для любой гладкой функции и и любого гладкого вектора У. 


р ди ди ы 
Здесь (Уи, И)! = = — — И. — кососимметричное скалярное произведение векторов Уи и И ,‚ причем 
х 


У ду 
(а,а)| =0 для любого вектора а. 
Применяя (24) к уравнению для плотности тока (4), получаем: 
"(Г = Еаг(П) + Ее"). (25) 


Используя формулу потоков (1), получим соотношение: 


_ Е ы 
Г(];) —-— ариг(С;Уф) — Б.У С;) + (СУМ), а (26) 
Отсюда с учетом свойств оператора Г: 
- Е - 2 
1: )= и — 5:0; (УС; Уф), + (УС, У) + С;т(У), =1,2. (27) 
- ОГ, 9, 
Учитывая (27) и г(Г) = = = Ди ‚ уравнение (25) можно записать: 
х 
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2 
ыы /Зщ |[ 2р А +12 ор, 0с, \0ф _ ар, А +2, 0, 00 |. 
ВТ ОХ @х / 0 ду ду ты 


дс д 0 Ос, ОУ, _ 0, 
+Е не ИЕ а— “+: У +Р (ас+5с,) | —^ 
Ох Ох ду . бх 0у 


Из (28) и (20) следует, что функция П определяется с точностью до константы, поэтому можно допу- 


(28) 





стить: 
1(й,0) =0. (29) 

Тогда из (23) получаем 

(А, =. (30) 

Условия (29) и (30) — краевые для функции 1. 

Математическая модель для гальваностатического режима в безразмерном виде. Для численного ис- 
следования краевых задач удобно перейти к безразмерному виду. Так можно упростить уравнения и выяснить 
фактическое количество и набор параметров, определяющих поведение системы. Безразмерные переменные 
описывают класс подобных процессов, характеризующихся одинаковым значением безразмерных чисел. 

Характерные величины, описывающие задачу. При моделировании процессов массопереноса в камере 
обессоливания ЭДА принимается ряд характерных значений: 

— для пространственных координат — межмембранное расстояние й; 


— для концентраций ионов — объемная концентрация электролита Со; 

— для скорости — средняя скорость вынужденного течения У, ; 

— для коэффициентов диффузии — коэффициент диффузии электролита р =р.О, ($, -5, )/ Ву, 
— для электрического потенциала — тепловой потенциал Фф, = КТ/ Иа 

— для плотности тока — величина &, = ЕРСу/й (аналог предельной плотности диффузионного тока); 
— для потока ионов — диффузионный поток ®=В9( ‚ й. 


Формулы перехода. Переведем уравнения в безразмерную форму с помощью следующих соотношений 
(индексом (и) обозначены безразмерные варианты величин): 


— 


ох У ЮГ У со 


Хх = у , 1 
й й й \ Со 
1 1 р к 
060 =. =, о р =. 
Фо й ' Л 
Система уравнений в безразмерной форме имеет вид (индекс (И) для упрощения записи опущен): 
1; =-=), С, Уф-Б, УС; + РеСУ, 1=12, (32) 
—ау =0, 1=Ь2, (33) 
дс 
Ат=-|| В, — р, 9% ыы гр ба :2р, 9% 2 |8 |+ 
де сх ду ты - Ох (35 
д д д д 0У, ду 
+Ре ча, — У, -Ре ры НМ У. +Ре(ас + 5,с,)| —^-—^ |, 
Ох Ох а. ду 0х ду 
ТЕ, + 5212 . (36) 
Система уравнений (29)—(35) содержит два безразмерных числа: Ре =Ууй/Рр — число Пекле и 


в й) ыы <“ 2 

&=8 „воВТ/ (Сий” Е”). Физический смысл параметра = состоит в том, что это удвоенный квадрат безразмерной 
ь | у, 

дебаевской длины р: &=2(1р /й)" [5]. Оценка величин параметров показывает, что при естественных для 


электродиализа условиях число Пекле имеет порядок 10? —10°, число & имеет порядок 10! —107” ‚то есть 


может считаться малым Параметром. 
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Для удобства численного решения преобразуем систему уравнений, подставив плотность потока (32) в 


уравнения (33) и (36): 
4’ (-=.0, С, Уф- Б,УС, +РеСУ)=0, #=1,2, (37) 
- 2 и 
ГУ: [ <), С. Уф-Б, УС, + Рес (38) 
= 


Таким образом, система уравнений содержит следующие неизвестные функции х, у: С, С2,Ф, [х,1.. 

Поля концентраций С\, С. и потенциала ф определяются решением уравнений (37), (34) соответственно. Ком- 

поненты плотности тока /,,[, вычисляются с помощью (38). Распределение скорости (5) в безразмерной фор- 
ме: 

У =0, У, =6х (1-х). (39) 


Присутствие малого параметра # в уравнении Пуассона (34) означает, что краевая задача является син- 
гулярно возмущенной. Это значительно усложняет ее численное решение, поскольку такие задачи относятся к 
жестким [26]. Потенциал электрического поля ф и концентрации ионов С, С> изменяются очень быстро в уз- 


ком пограничном слое, толщина которого равна длине Дебая [> [5]. Для решения этой проблемы целесообраз- 
› р 


но уплотнить вычислительную сетку в пограничном слое и использовать специальные методы решения жест- 
ких задач [26]. 
Граничные условия в безразмерной форме. На межфазной границе «АОМ — раствор» (х = 0): 





ас! дФ (-7>)1,(0, у) 
И + С ЕЕ 0, р ЗЕ ь 40 
Ох т Ох } У) 52 ) 
С> (0, у)= М. (41) 
Ф(0,у)= 0, (42) 
21 (0,5) =0 (43) 
ОХ 
На границе «раствор — КОМ» (х = |): 
С(Ь У) = №, (44) 
9С> 0Ф (7: )7,,(, у) 
А 45 
Ох 622 ОХ } У) 5205 ) 
— О" +; р ОСТ +2 р ОС 
дФ ду ах “^^ а 
(= о (, >) (46) 
2 2 
Ох &1 ВАС1 + 5262С> 
21, у)=0. (47) 
ОХ 
На входе в канал (у = 0): 
С: (х0)=1 #=12, (48) 
МТ [рр ры 49) 
ду = +5505 ду ду 
1(х.0)= 0. (50) 
На выходе из канала (у = [): 
- р 5;С; 5 ый = 0, = 1,2, (51) 
ду ду 
м = 0, (52) 
ду 
1(х,1) =-, у = с0п5 . (53) 


После численного расчета системы (34), (35), (37)—(53) скачок потенциала Аф в канале обессоливания 


определяется по формуле: 
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й 17 
Ар = | ®(ьу)4у. (54) 
0 


Численное решение найдено методом конечных элементов с помощью пакета Соп1$01 Ми!арВу$1с$ на 
неравномерной вычислительной сетке (плотность элементов сетки увеличена у границ «раствор — мембра- 
на») [27]. 


Результаты исследования. Вычисления проведены для &=1,9.10”, Ре= 2355 , что соответствует 


следующим значениям параметров системы: 
— входная концентрация раствора электролита МаС Си = 0,1 моль/М’; 
— температура 7 `= 298 К; 
— коэффициенты диффузии катионов и анионов соответственно О; = 1,33:10`° м”/с, 2» = 2,05:10° м”/с; 
— числа переноса противоионов в мембранах Т! = 0,972, Т. = 1; 
— зарядовые числа ионов 1 = 1, > =-1[; 
— отношение концентрации противоионов на границе с мембранами к ее значению на входе в канал № = М, = 1; 
— й= 10 м — ширина канала: 
— 1=2.10 м — длина канала; 
— И = 3, 8-10 м/с — скорость прокачки раствора. 
На рис. 2 приведены поля концентраций С! и С›, потенциала ф и функции 1, рассчитанные при плотно- 
сти тока („= 1,5йт» где пм — это предельная плотность тока, определенная по формуле Левека (в безразмерной 
форме) [28]: 


1/3 


2 
. 24 =0,2. |. (55) 


т — › 7 
1 -& 





Здесь й = 0,395 — число переноса катионов в растворе [9]. 
й 








в) г) 


Рис. 2. Поля концентрации катионов С! (а) и анионов С. (6), потенциала ф (в) и функции т (г). 


Расчет по модели для гальваностатического режима при плотности тока й, =Ъ5й„ 
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В канале действует вынужденное течение, поэтому в областях у мембран истощение концентрации 
ионов увеличивается при удалении от входа в канал (вдоль направления, тангенциального к поверхностям мем- 
бран). Соответственно, толщина расширенной области пространственного заряда увеличивается вдоль кана- 








ла (рис. 3). 
01 
$ \2 1 
® 
5 0,05 
С. 
ыы 
х 0,98 0,99 1,00 
х 
а) 6) 
0,01 
= 0,005 
О 
0,96 0,98 100 


ви 
в) 
Рис. 3. Результаты численного расчета по модели для гальваностатического (сплошные линии) и потенциостатического 


(пунктирные линии) режимов при плотности тока Г, =1,51„, в сечениях у = 0, (1), у=0,41 (2), у=0,91(3): 


концентрационные профили С! и С. (а); увеличение рис. 3, а (6); плотность пространственного заряда р = <1С] + <2С. (в) 


Вольтамперная характеристика (ВАХ, кривая | на рис. 4) рассчитана на основе гальваностатической 


модели. При этом задаваемая плотность тока изменялась: („= 0; 0,0051; ...; 1,ЭЙрт. 


40 





30 
= 
5 
— = 
д к. 
— 20 
э 
= 
В 
о 10 
> 
3 
с 
> 
== 


О 10 20 30 40 
Ао 
Рис. 4. ВАХ, рассчитанные по моделям для гальваностатического (1) и потенциостатического (2) режимов; 
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Предельная плотность тока определена по точке пересечения касательных к ВАХ в начальной части и 
на наклонном плато [5]. Здесь она совпадает со значением й„, оцененным по приближенной формуле Леве- 
ка (55). Также ВАХ была рассчитана на основе потенциостатической модели при изменении скачка потенциала 
Дф = 0; 0,4; ...; 40 (кривая 2 на рис. 4). Как видно из рис. 4, ВАХ Ти 2 совпадают. Таким образом, существует 
однозначное соответствие каждому скачку потенциала некоторой плотности тока, и наоборот. 

Рис. 3 и4 демонстрируют достаточно хорошее совпадение различных физико-химических характери- 
стик переноса, рассчитанных по моделям для гальваностатического и потенциостатических режимов. Это дока- 
зывает адекватность построенной авторами модели переноса в гальваностатическом режиме. 

Вычислительные затраты предлагаемой модели оценивались следующим образом. Фиксировалось вре- 
мя, затраченное на решение краевых задач для гальваностатического и потенциостатического режимов. При 


этом значения задаваемой точности вычислений, параметры системы были одинаковыми, а 1, = 1,5 Г„. Затем 


показатели времени сравнивались. Таким образом выяснилось, что расчет по гальваностатической модели тре- 
бует в 1,6 раз больше времени. Это связано с тем, что: 
— гальваностатическая модель содержит дополнительное уравнение для определения распределения плотности 
тока; 
— для потенциала на границе х = 1 установлено условие второго рода (46). 

Комплексный потенциал электромембранной системы для гальваностатического режима. В рабо- 
те [25] предложено обобщение импеданса электрохимической системы с помощью функции \, введенной в 


условиях электронейтральности. Аналогичное обобщение приемлемо и в данном случае. Отметим, что при ис- 
пользовании электрохимического импеданса объект рассматривается только как «черный ящик», и его внут- 
ренние свойства определяются косвенно. Полученные выше результаты позволяют ввести понятие комплексно- 
го потенциала электромембранной системы: Р=Фф-+1.1. Комплексный потенциал — это функция координаты 


любой точки внутри объекта, поэтому, в отличие от электрохимического импеданса, он позволяет исследовать 
внутренние свойства объекта. 

Обсуждение и заключения. Описан новый метод математического моделирования стационарного 
процесса массопереноса в гальваностатическом режиме для мембранных систем. При этом рассмотрен двумер- 
ный случай с использованием специального граничного условия, позволяющего задать плотность тока в систе- 
ме. Представлены уравнения для электрической функции тока. Результаты численного решения по потенцио- 
статической и гальваностатической моделям хорошо согласуются. Это показывает адекватность предлагаемой 
модели переноса в гальваностатическом режиме. 

Разработанная модель позволяет интерпретировать результаты экспериментальных исследований пере- 
носа ионов в мембранных системах, если данный процесс протекает в гальваностатическом режиме. Некоторые 
электрокинетические процессы связаны с появлением расширенной области пространственного заряда при 
сверхпредельных токах. Описывая формирование указанной области, можно выяснить, каким образом завися- 
щие от нее процессы влияют на перенос ионов в гальваностатическом режиме. 
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